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Предложена принципиально новая кинетическая модель процесса дегидрирования н-парафинов С 10 ~Св, в основе которой ле- 
жит формализованный механизм превращения углеводородов на поверхности РІ-контакта. Модель учитывает дезактивацию 
катализатора вследствие коксонакопления с учетом специфики перерабатываемого сырья и условий ведения процесса. 


В мире из общего потребления нормальных па- 
рафинов 75 % используется для производства ли- 
нейных алкилбензолов (ЛАБ) - сырья для синтети- 
ческих моющих средств. Потребность в синтетиче- 
ских моющих средствах будет расти и в третьем ты- 
сячелетии. Поэтому повышение эффективности 
процесса производства ЛАБ, в частности, дегидри- 
рования парафинов на Рі-катализаторах как нео- 
тъемлемой его части, является важной и актуаль- 
ной задачей [1]. Однако существует ряд факторов, 
усложняющих ее решение. В первую очередь это 
многостадийность и многокомпонентность систе- 
мы. Во-вторых, это нестационарность процесса, 
протекающего в условиях постепенного снижения 
активности катализатора в результате накопления 
на его поверхности коксовых отложений. Решение 
этой проблемы может быть реализовано с исполь- 
зованием метода математического моделирования. 


Исследования последних лет показывают, что с 
помощью технологических приемов и придания 
определенных свойств катализаторам в ряде случа- 
ев удалось существенно усовершенствовать промы- 
шленные каталитические процессы [2]. 

Так, основными методами повышения эффек- 
тивности процесса получения н-моноолефинов до 
настоящего времени остается разработка и внедре- 
ние новых катализаторов дегидрирования н-парафи- 
нов. При этом решаются, как правило, две задачи: 
создание, промышленные испытания и выбор ак- 
тивного, высокопроизводительного катализатора; 
повышение качества продукта за счет увеличения се- 
лективности и стабильности выбранного контакта. 

Все перечисленные выше приемы относятся к 
типу экспериментальных (инструментальных) и 
требуют, в свою очередь, крупных экономических 
вложений и длительного времени. 


81 


Известия Томского политехнического университета. 2006. Т. 309. № 5 


В настоящее время все большее распростране- 
ние находят методы математического моделирова- 
ния, реализованные в виде технологических моде- 
лирующих систем, которые позволяют прогнози- 
ровать длительность эксплуатации и срока регене- 
рации используемого катализатора, оптимизиро- 
вать технологический режим с учетом динамики 
процесса и рассчитывать ряд других важных пока- 
зателей промышленной технологии. В основе раз- 
рабатываемых моделей должны лежать физико-хи- 
мические закономерности превращения углеводо- 
родов на поверхности катализатора. 

Таким образом, целью данной работы является 
разработка формализованного механизма превра- 
щения веществ в процессе получения н-монооле- 
финов с учетом дезактивации катализатора вслед- 
ствие его закоксовывания. 

Кинетические модели разработаны достаточно 
полно лишь для смесей из нескольких компонен- 
тов. В этом случае не возникает принципиальных 
затруднений при построении уравнений математи- 
ческой модели, а также с определением параметров 
путем решения обратной кинетической задачи. 
Для многокомпонентных промышленных процес- 
сов, одним из которых является процесс дегидри- 
рования парафинов С 10 — С 13 данные методы непри- 
годны, т. к. кроме многокомпонентности в этих си- 
стемах имеет место и многостадийность реакций. 
При этом отсутствие информации о концентрации 
промежуточных компонентов усложняет задачу 
построения адекватной реальному процессу неста- 
ционарной кинетической модели. 

В связи с этим, следуя по пути агрегирования 
компонентов для описания сложной кинетики 
превращения веществ в ходе процесса получения 
н-моноолефинов С 10 — С 13 , необходима формализа- 
ция механизма при обосновании уровня детализа- 
ции, достаточном для описания кинетики сложных 
углеводородных смесей с учетом специфики при- 
меняемых катализаторов. 

Крайним случаем формализации механизма яв- 
ляется принцип «черного ящика», когда известны 
лишь входные и выходные параметры системы, а 
их взаимосвязь определяется некоторой статисти- 
ческой функцией. Однако в этом случае не заложе- 
ны физико-химические закономерности превра- 
щения веществ в ходе процесса. Кинетические па- 
раметры определяются почти вслепую (отсюда и 
название «черный ящик») [3]. Тем не менее, прин- 
цип «черного ящика» широко используется для 
анализа сложных систем, которым является любое 
химическое производство. Математические моде- 
ли, построенные на основе этого принципа, не ин- 
вариантны к изменению качества сырья и техноло- 
гических условий получения продуктов, т. к. в ос- 
нову этих моделей не положены физико-химиче- 
ские закономерности превращения углеводородов. 
На практике модели, которые построены на основе 
принципа «черного ящика», позволяют выработать 
некоторую рабочую гипотезу или предположение о 


существовании закономерностей превращения 
сырья в продукт. 

Вторым предельным способом представления 
информации является составление подробной систе- 
мы показателей, которые функционально связыва- 
ют закономерности превращения сырья в продукт. 
Ввиду сложности данной системы, трудностью, а 
иногда и невозможностью определения ее параме- 
тров, получить конечный результат затруднительно. 

Разработанный и развиваемый нами подход фор- 
мализации механизма превращения веществ осно- 
ван на физико-химических и технологических зако- 
номерностях протекания реакций в промышленных 
условиях. При использовании данного подхода из 
рассмотрения исключаются те закономерности, ко- 
торые не оказывают существенного влияния на тече- 
ние процесса, а также те закономерности, численное 
значение параметров которых трудно определить. 
При использовании формализованного способа есть 
возможность при некоторых допущениях изучить 
аппарат и явления, протекающие в нем. 

Описание механизма сложных многокомпо- 
нентных реакций на Рі-катализаторах - начальный 
этап компьютерного прогнозирования и постро- 
ения физико-химической модели. При этом осно- 
ву составляет термодинамика процесса. Результа- 
том термодинамического анализа является список 
реакций, которые, вероятно, протекают в заданных 
условиях. Для количественного описания скоро- 
стей протекающих реакций необходимо формали- 
зовать - выбрать список основных компонентов и 
реакций. С физико-химической точки зрения 
объединению должны подлежать вещества и реак- 
ции с подобными, с точки зрения кинетики, свой- 
ствами. С математической точки зрения формали- 
зованный механизм должен включать реакции, для 
которых возможна количественная оценка кинети- 
ческих и термодинамических параметров [4]. 

Так, реакция дегидрирования высших н-пара- 
финов протекает на металлическом центре моди- 
фицированных алюмоплатиновых катализаторах 
под небольшим давлением Н, при температуре 
450. ..500 °С. В этих условиях выходы олефинов 
ограничены термодинамическим равновесием и 
составляют, как правило, 7. ..11 % [5]. 

Образование олефинов является первой стади- 
ей ароматизации парафинов [6]. Эта реакция на 
оксидных катализаторах протекает в основном по 
консекутивному механизму по схеме парафин—» 
олефины—» диолефины— »триолефины-»циклогек- 
садиены— »ароматические углеводороды [7]. В де- 
гидрогенизатах высших парафинов обычно об- 
наруживаются кроме олефинов диены и аромати- 
ческие углеводороды, а также продукты крекинга и 
скелетной изомеризации. Если отщепление водо- 
рода и образование олефинов, диенов и триенов 
происходит на дегидрирующих металлических 
центрах катализатора, то крекинг и изомеризация, 
по-видимому, идут на носителе; циклизация три- 
енов может протекать термически. 
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Таблица. Брутто-реакции процесса дегидрирования н-парафинов 


Дегидрирование парафинов до 
моноолефинов 


Дегидрирование моноолефинов до 
диолефинов 


с 1 П н„ 


с„н 


±С 1П Н, П +Н, 


11 Г1 24 <- 


с„н 


±с„н,,+н. 


12^26 


С„Н 


:±С 12 Н 24 +Н 2 


13 Г1 28 <- 


— С 13 Н 26 +Н 2 


с,„н. 


с„н. 


с„н 


^С 1П Н 10 +Н, 


^С„н, п +н, 


12 Г1 24 


±С, 2 Н„+Н 2 


с„н 


13 П 2б <- 


— С 13 Н 24 +Н 2 


Изомеризация 


Дегидроциклизация 


п-С 10 Н 22 ^=^і-С 10 Н 22 

«12 и 

п 'СцН 24 <п==±і-С 11 Н 24 


П ‘С 12 Н 26 <===^і-С 12 Н 26 

кто,, 


п-С 13 Н 28 ^=±і-С ІЗ Н 28 


К 


п-С 10 Н 18 ^±С 6 Н 5 -С 4 Н 9+ 2Н 


п-С„Н 




Ь±С 6 Н 5 -С 5 Н п + 2Н, 


п-С 12 Н 22 С 6 Н 5 - С 6 Н 13 + 2Н 2 

к-ыц 

п-С 13 Н 24 ^==> С 6 Н 5 - С 7 Н 15 + 2Н 


Гидрокрекинг парафинов 


с, п н„+н 


1 0 П 22 т П 2 4- 


— С 5 Н 12 +С 5 Н 12 


СцН 24+ Н 2 ф=2=±С 5 Н 12 +С 6 Н 14 
^■12^26 Н 2 4 
СізН 28 + Н 2 4 


-С 6 Н 14 +С 6 Н 14 


— С 6 Н 14 +С 7 Н 16 


На алюмоплатиновых катализаторах процесс аро- 
матизации парафинов протекает более сложно, так 
как в их присутствии могут образоваться пятичлен- 
ные цикланы [8]. Однако наличие в катализаторе де- 
гидрирования специальных добавок подавляет эту 
реакцию, и процесс ароматизации парафинов в усло- 
виях дегидрирования на Рі-катализаторе осущест- 
вляется по тому же механизму, что и на оксидных ка- 
тализаторах [9-1 1]. При дегидрировании высших н- 
парафинов на всех катализаторах получаются олефи- 
ны в основном с внутренними двойными связями. 

Количественно оценить энтропийные и энталь- 
пийные характеристики промежуточных стадий за- 
труднительно, т. к. невозможно измерить концен- 
трацию веществ, участвующих в этих стадиях. На- 
бор брутто-реакций процесса дегидрирования н- 
парафинов в формализованном механизме пред- 
ставлен в таблице, где к, - константы скоростей хи- 
мических реакций. 

Таким образом, с учетом формализации основ- 
ных стадий процесса дегидрирования н-парафинов 
на Рі-катализаторе была выбрана следующая схема 
превращений веществ, рисунок. 


Основными реакциями, протекающими в ходе 
данного процесса, являются: 

• дегидрирование н-парафинов до н-олефинов 
(целевая реакция); 

• образование диолефинов из олефинов; 

• образование ароматических углеводородов из 
образовавшихся диолефинов; 

• расщепление (крекинг) углеводородов с образо- 
ванием продуктов уплотнения и легких углево- 
дородов; 

• реакция коксообразования; 

• изомеризация. 

На основании приведенной схемы превраще- 
ния веществ была разработана кинетическая мо- 
дель процесса, которая позволила количественно 
оценить скорости превращения углеводородов в 
зависимости от времени реакции, состава исход- 
ной смеси, температуры и давления в системе. Учет 
макрокинетики позволил нам рассмотреть химиче- 
ские стадии на поверхности Рі-катализатора в их 
взаимосвязи с процессами энерго- и массоперено- 
са в промышленном аппарате определенной кон- 
струкции, т. е. с учетом гидродинамической соста- 
вляющей в модели реактора. 
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Рисунок. Схема превращения углеводородов процесса получения н-олефинов 


Так как промышленные процессы получения 
н-моноолефинов осуществляются в адиабатиче- 
ских реакторах идеального вытеснения, где длина 
реактора много больше его диаметра (Ь»В), то из- 
менение концентрации /-го компонента в таком 
реакторе по времени запишется уравнением: 


дС , 
ді 


-ІІ^ + О-еХ, 
ді 


( 1 ) 


где С ( - концентрация /-го углеводорода, моль/м 3 ; 
II - линейная скорость потока, м/ч; г, - скорость 
реакции, моль/м 3 -ч; е - порозность слоя катализа- 
тора; / - время, ч. 

Время пребывания реагентов в реакционной зоне: 

х=ьт=ѵіс, 


где С - часовой расход сырья, м 3 /ч; V- объем ката- 
лизатора, м 3 . 

В этом случае уравнение (1) перепишется следу- 
ющим образом: 


дС < г^ дС , п ^ 

— +С — - = (1 ~е)г 
ді дѴ ' 


В условиях нестабильной нагрузки установки 
по сырью используется «приведенное время» или 
суммарный объем переработанного сырья после 
регенерации катализатора с=Сі, м 3 . 

Тогда 


дг дѴ 




Начальные и граничные условия: 
г = 0 : С, = 0, 

Ѵ = 0 :С, = С„. 


Уравнение теплового баланса в дифферен- 
циальной форме запишется как: 


У А Н.г, 

дТ дТ 

— + С — = -(1-е)- 

ді дѴ с р 


Здесь Т- температура процесса, К; АД - тепло- 
вой эффект реакции, Дж/моль; с р - теплоемкость 
смеси, Дж/моль-К; р - плотность смеси, кг/м 3 . 

Или в приведенной форме: 



дѴ 


Ха н, Гі 
-(і -<0- 

с ,>р 


Начальные и граничные условия: 

г=0: Т= 0, 

К=0: Т=Т вх . 

При составлении данной системы приняты сле- 
дующие допущения: 

• Отсутствие диффузионных осложнений. 

• Порядок реакции по /-му компоненту соответ- 
ствует его стехиометрическому коэффициенту. 

• В качестве гидродинамического режима выбран 
режим идеального вытеснения. 

Нестационарный характер протекания промы- 
шленного процесса обусловлен изменением актив- 
ности катализатора в ходе его эксплуатации. Отло- 
жение продуктов уплотнения, блокирующих ак- 
тивные центры катализатора, является существен- 
ной причиной падания его активности. 

Кокс, оседающий на Рі-контакте при дегидри- 
ровании индивидуальных углеводородов, тормозит 
основную реакцию. Существует определенная связь 
между каталитической активностью, количеством 
кокса и временем, в течение которого катализатор 
находился в контакте с углеводородом. Оседающий 
кокс блокирует часть поверхности, вызывая диффу- 
зионное торможение основной реакции. Однако 
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имеются различия во взглядах, как блокируется по- 
верхность катализатора. В ряде работ описаны под- 
ходы к моделированию каталитических процессов, 
учитывающие дезактивацию катализатора из-за 
коксоотложения [5, 12, 13]. Большинство извест- 
ных уравнений кинетики образования кокса пре- 
дельно формализованы, они содержат большое чи- 
сло эмпирических коэффциентов. Математические 
модели дезактивации катализаторов в основном 
статистические, настроены на узкий интервал из- 
менения параметров и практически не учитывают 
концентраций реагирующих веществ. 

Изменение каталитической активности проис- 
ходит за счет изменения числа активных центров 
катализатора. Для количественного определения 
активности катализатора вводится параметр а(і), 
который соотносит концентрацию активных цен- 
тров катализатора в любое время с начальной кон- 
центрацией центров: 


а(і) 


т 

щУ 


( 2 ) 


где 5(?) — количество центров на закоксованном 
катализаторе; У(? 0 ) - общее количество центров на 
свежем катализаторе. 

В качестве допущений можно принять, что кокс 
образуется преимущественно за счет поликонден- 
сации ненасыщенных моноциклических углеводо- 
родов. Предшественники кокса получаются из 
углеводородов, находящихся в газовой фазе. Таким 
образом, число активных центров на катализаторе 
зависит от количества кокса, в который превраща- 
ются углеводороды из газовой фазы: 


г/5(7) 

йС к 




V і= 1 


К,(х)С,(х) 5(() 


( 3 ) 


или 
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V *=і ) 


а )■ 


(4) 


Здесь К,(х) - коэффициент скорости блокиров- 
ки активных центров у-ой реакции из С : (х) компо- 
нентов; С{х) - концентрация компонента /-ой го- 
могенной группы, содержащего х атомов углерода в 
углеводородной цепи; С к - концентрация кокса на 
катализаторе. 


Интегрирование ур. (4) дает зависимость изме- 
нения активности катализатора в данном типе реак- 
ций от общего содержания кокса на катализаторе: 

а ] =Ае а,Ск , (5) 

где а, - коэффициент отравления, зависящий от 
соотношения кислотного и основного фактора; А - 
энтропийный фактор процесса. 

Для определения численных значений ис- 
пользуются результаты лабораторных испытаний 
или промышленных испытаний катализатора. 

Скорость у'-ой химической реакции /-го компо- 
нента с учетом дезактивации катализатора предста- 
вляется через а/. 

П -С.к а,. (6) 

Здесь ^ - константа скорости реакции на свежем 
катализаторе. Различие скорости закоксовывания 
разных центров катализатора в одном и том же про- 
цессе отражено в значениях а г Процесс коксообра- 
зования затухает с течением времени, что объясня- 
ется уменьшением числа свободных центров. Это 
соответствует математическим выражениям (2) — (5). 

Уровень активности а і можно рассматривать как 
изменение удельной поверхности катализатора в 
каждый момент времени. Ввиду трудности опреде- 
ления числа центров эффективен способ их кос- 
венного учета через адсорбцию (6). 

Изменение технологических параметров влияет 
на отравление активных центров катализатора кок- 
сом. Изменение температуры, давления, кратности 
циркуляции водородсодержащего газа могут приве- 
сти как к равномерной дезактивации, так и к «пере- 
распределению» кокса между активными центрами. 

Таким образом, разработан формализованный 
механизм превращения углеводородов на Рі-ката- 
лизаторах дегидрирования н-парафинов и предло- 
жена схема превращения углеводородов в ходе дан- 
ного процесса. 

Построена нестационарная кинетическая мо- 
дель промышленного процесса получения н-оле- 
финов с учетом закономерностей изменения ак- 
тивности катализатора в зависимости от химиче- 
ского состава сырья и условий ведения процесса, а 
также с учетом дезактивации катализатора вслед- 
ствие его закоксовывания. 
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